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Low-field (6 — 110 mT) magnetic resonance of bismuth (Bi) donors in silicon has been observed by 
monitoring the change in photoconductivity induced by spin dependent recombination. The spectra 
at various resonance frequencies show signal intensity distributions drastically different from that 
observed in conventional electron paramagnetic resonance, attributed to different recombination 
rates for the forty possible combinations of spin states of a pair of a Bi donor and a paramagnetic 
recombination center. An excellent tunability of Bi excitation energy for the future coupling with 
superconducting flux qubits at low fields has been demonstrated. 



Among a variety of qubits investigated for the realiza- 
tion of solid-state quantum computers, superconducting 
qubits are the leading candidates for quantum proces- 
sors because of their fast operation capabilities (tv/2 pulse 
shorter than 10 ns)ii However, the shortcoming of their 
relatively fast decoherence time needs to be overcome 
by connecting to memory qubits that can store quan- 
tum information throughout the course of computation. 
This requires memory qubits working under low magnetic 
field, typically below 10 mT for aluminum superconduct- 
ing qubits? 2 - since they cannot operate at higher fields. 

In this context, the bismuth (Bi) donor in silicon (Si) 
has attracted much attention recently. Its large hypcrfinc 
interaction a/h = 1.4754 GHz (Ref. Q) and the 209 Bi 
nuclear spin 1 = 9/2 give a large zero- field splitting of 7.4 
GHz. This splitting is comparable to the typical energy 
splitting between \R) and \L) states of superconducting 
flux qubitsii^ Thus, coupling between a Bi spin qubit 
and a superconducting flux qubit on Si is in principle 
possible via a microwave photon through a waveguide. 
A proposal of such an application 5 has prompted exten- 
sive research on the Bi donor spins in Si very recently££ 
Starting from the spectroscopic analysis of the electron 
paramagnetic resonance (EPR),»i the electron spin re- 
laxation time TijiliS decoherence time T2r^& and super- 
hyperfine interaction with nearby 29 Si nuclear spins^ were 
investigated. Moreover, the coherent transfer between 
electron and 209 Bi nuclear spins 6 and dynamic nuclear 
polarization of 209 Bi£i£ were achieved. Yet all of these 
EPR studies were performed at 9 GHz (around 320 mT) 
and at 240 GHz (around 8.6 T) excitation frequency. 

In this paper we report on low-field (6 — 110 mT) radio 
frequency (20 - 400 MHz) and microwave (8.141 GHz) 
magnetic resonance, as well as X-band (9 GHz) magnetic 
resonance, of ion-implanted Bi donors in Si using a highly 
sensitive, spin dependent recombination based magnetic 
resonance (SDR-MR) method^-i^ 



The samples were prepared from highly resistive (> 
3000 f2-cm), float-zone (FZ) n-type silicon wafers im- 
planted with Bi ions at room temperature with a total 
fluence of 2 x 10 12 cm -2 . The implantation energies are 
300 and 550 keV with the doses of 0.7 x 10 12 and 1.3 x 10 12 
cm" 2 , respectively. This condition leads to the Bi con- 
centration of about 1.8 x 10 17 cm -3 in the depth of 90 to 
150 nm from the surface. The post-implantation anneal- 
ing was performed at 650 °C for 30 min in an evacuated 
quartz tube. This annealing condition has been shown 
to achieve the electrical activation of 80 % ) Z i 14 i 15 thus 
yielding about 4.8 x 10 11 Bi donors in a 0.3 cm 2 sample 
area. The SDR-MR spectra were recorded at 16 K with 
a commercial continuous wave EPR spectrometer (JEOL 
JES-RE3X) working at X-band (9 GHz microwave) with 
a homemade coil for radio frequency (20 — 400 MHz) and 
microwave (8.141 GHz) irradiation to induce magnetic 
resonance at low field (6 — 110 mT). Continuous illumina- 
tion with a 100-W halogen lamp generates photocarriers 
in the sample. Magnetic resonance can enhance the spin- 
dependent recombination, which decreases the density of 
photocarriers. Then, the absorption of the microwave 
electric field by the photocarriers is decreased, leading 
to an enhancement in the Q-factor of the cavity. Thus, 
the effect of magnetic resonance can be detected simply 
as the change in the X-band microwave reflection from 
the cavity. The second derivative of the reflected inten- 
sity with respect to the field modulation was recorded as 
an SDR signal to reduce the broad cyclotron resonance 
lines and the background change of the sample resistiv- 
ity during the magnetic field scan. Note that because 
of our high power (80 mW) saturating excitation, i. e., 
making the populations of the ground- and excited-states 
the same, the conventional EPR absorption signal is sup- 
pressed. Such saturation is necessary to flip one of the 
spins in a pair of Bi and defect to induce SDR as we will 
discuss later. 
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The spin system of an isolated Bi donor in static mag- 
netic field B (electron spin S = 1/2 and 209 Bi nuclear 
spin / = 9/2) can be represented by the spin Hamilto- 
nian: 



H = g^BB Q S z - g n ^ n B I z +aS ■ I, 



(1) 



where fis and fi n are the Bohr and nuclear magnetons, 
and g = 2.0003 (Ref. 0) and g n = 0.914 (Ref. @) are 
Bi electron and nuclear ^-factors, respectively. The SDR 
method requires a coupled pair of electron spins In 
the present study, the partner of the Bi donor electron 
spin is the electron spin of a deep paramagnetic recom- 
bination center (PRC) which is supposed to be created 
during the implantation process and not completely re- 
moved by controlling the annealing conditions^ We have 
attempted to identify the symmetry of the deep PRC by 
tracing the angular dependence of the EPR peaks. How- 
ever, the peaks were too broad to draw conclusions. 

The X-band (9.076 GHz) SDR-MR spectrum measured 
without radio frequency excitation is presented in Fig. 
1(b). The peaks labeled as Bi-X and PRC indicate ten 
EPR- "allowed" transitions of the Bi donors and one EPR 
transition of a PRC, respectively, corresponding to the 
intersections with the 9.076 GHz excitation in Fig. 1(a). 
Here the EPR transition frequencies of the Bi donor 
(solid curves )22 were calculated as functions of the static 
magnetic field Bq using the spin Hamiltonian in Eq. (1) 
and that of PRC (dashed curve) was calculated as an 
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FIG. 1. (a) Calculated EPR transition frequencies of Bi 
donors (solid curves) and deep paramagnetic recombination 
center (PRC) (dashed line) with 9.076 GHz microwave exci- 
tation frequency (dot-dashed line). The intersections of the 
PRC transition energy with bismuth donor transitions (open 
triangles) and with the 9.076 GHz microwave (open square) 
are also shown, (b) An SDR-MR spectrum of Bi donors in 
Si recorded at 16 K under illumination. The 9.076 GHz mi- 
crowave is used both to induce Bi EPR transition and to probe 
the change in the sample photoconductivity. CR1 and CR2 
are the cross relaxation signals between Bi and PRC. The 
open circles in (b) indicate simulated intensities using the 
SDR model described in the text with the parameter value of 
flp/i?ap = 0.01. 



isotropic, nuclear spin free, paramagnetic center S = 1/2 
and g » 2.005(3). The same notation as in Ref. [20| 
for labeling Bi eigenstates is used; the labels 1 to 20 in 
increasing order of energy. 

In addition, two lines labeled as CR1 (186 mT) and 
CR2 (259 mT) arise due not to the resonance with the 
9.076 GHz microwave but to cross relaxation (CR) be- 
tween particular Bi donor transitions and the PRC tran- 
sitions, in a way very similar to the cross relaxation be- 
tween phosphorus donors and SL1 centers in Si observed 
by electrically detected magnetic resonance.— This as- 
signment of CR1 and CR2 is further justified in Fig. 
2. Even with the decrease in the microwave excita- 
tion frequency, the position of these lines remains the 
same whereas the Bi EPR line positions shift to lower 
fields. This proves the presence of coupling between the 
Bi donor and PRC electron spins, which is requisite for 
the SDR detection method. 

We should emphasize that even at the same X-band 
resonance of the Bi donor, the observed SDR-MR line in- 
tensity distribution is clearly different from that observed 
in the conventional EPR spectra.^^^ The intensity dif- 
fers for the ten different transitions in the present SDR- 
MR whereas it is practically the same in conventional 
EPR, reflecting simply the thermal equilibrium popula- 
tion difference between the involved levels. Furthermore, 
the observed line-dependence of the SDR-MR intensity 
is distinctively stronger than the line dependence in the 
EPR transition^ 

Figure 3(a) shows the SDR-MR spectra probed by the 
same X-band microwave but with additional radio fre- 
quency excitation of 50 or 200 MHz. A simulation of the 
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FIG. 2. Cross relaxation lines CR1 (a) and CR2 (b) together 
with SDR-MR Bi lines detected at different microwave fre- 
quencies. The Bi EPR lines labeled as 1 16) -O- |5) in (a) and 
1 14} -R- 1 7) in (b) shift with the resonant frequency whereas 
the lines CR1 and CR2 do not. Additionally, the CR1 line 
is narrower than CR2 because, as shown in Fig. 1(a), the 
difference in the field-derivative of the transition frequency 
between the paramagnetic recombination center and resonant 
Bi transition is larger. 
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FIG. 3. (a) Low-field SDR-MR spectra with 50 and 200 
MHz resonance frequencies together with simulation of the 
200 MHz spectrum. The transitions by the radio frequencies 
and X-band 9.076 GHZ microwave are labeled as Bi-RF and 
Bi-X, respectively. The line intensities are simulated using 
the same model and parameters as for Fig. 1(b). (b) The 
Bi-RF line positions observed at various resonant frequencies 
(solid circles) together with calculated resonant fields (solid 
lines). 



200 MHz SDR-MR spectrum based on the SDR model 
that will be introduced later is also shown. Figure 3(b) 
shows the observed SDR-MR line positions (solid circles) 
for radio frequency excitation ranging from 20 to 400 
MHz together with the calculated magnetic resonance 
transition frequencies for the Bi donor (solid curves). All 
of these simulated Bi magnetic resonance transitions are 
the 209 Bi NMR transitions in the high field limit. While 
the number of calculated Bi NMR lines in Fig. 3(b) ap- 
pears ten, all but the lowest- and highest-field lines are 
nearly-degenerate doublets, i.e., each is composed of two 
lines separated by exactly twice the nuclear Zeeman split- 
ting energy. Therefore, among the eighteen EPR transi- 
tions of the Bi donor excited by these radio frequencies, 
we observed clearly the two non-degenerate lines labeled 
as Bi-RF in Fig. 3. The remaining eight doublets are 
too weak to be observed with the current experimental 
conditions. 

As mentioned above, considering only the EPR transi- 
tion probabilities^ cannot describe the intensity of SDR- 
MR. In the following, we shall present a model, based on 
the SDR model developed in Refs. [22T - [24j ] . to simulate 
the SDR-MR spectra that are shown in Fig. 1(b) and Fig. 
3(a). The SDR signal intensity, measured by probing the 
microwave intensity reflected by the cavity, is linear to 
the sample photoconductivity a. Its change by magnetic 



Act oc - 2J Ri,n 



[N itll {w^oo)-N i!fi (w = 0)}, (2) 



where the subscripts i and /i denote the Bi donor and 
PRC spin states, respectively. Ri tli and 7V ijM are the re- 
combination rate and population of the specified pair. 
Here the square bracket represents the change in the pop- 
ulation from off-resonance (w = 0) to saturated magnetic 
resonance (w — > oo) conditions where w is the excitation 
power. Eq. (2) is valid when the recombination rates can 
be assumed dominant over the pair generation and dis- 
sociation rates as well as the spin-lattice relaxation and 
spin decoherence rates. Furthermore, if only one transi- 
tion between two Bi states i and j is selectively excited, 
the change in photoconductivity becomes: 

Aa(i,j) w -5^[i2*, M JV ilM (oo) + Rj^N j>IM (oo) 

-(R^N^(0) + R^N^(0))}. (3) 

Using the rate equations described in Ref. [17], Eq. (3) 
simplifies to: 
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(4) 



The first term of Eq. (4) corresponds to the number 
of recombining pairs when the resonance is saturated 
whereas the second and third terms are the off-resonance 
terms. Thus, a large change in SDR signal should be ob- 
tained when either Ri^ or Rj t n is much smaller than the 
other. Finally, to evaluate the recombination rates Ri^ 
and Rj^ , the product state of the Bi donor and the PRC 
is considered: 



|*> |a*> =(cos4(f? )|l/2,M-l/2) 

+ sm<f> i {B ) |-l/2, M+ 1/2)) \fx) . 



(5) 
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FIG. 4. Low-field SDR-MR spectra probed by the 9.076 GHz 
X-band reflection (a) with a single 8.141 GHz excitation fre- 
quency and (b) with the same 8.141 GHz excitation plus an 
additional 100 MHz radio frequency. The lines resonant with 
the 8.141 GHz microwave and the 100 MHz radio frequency 
are labeled as Bi-8 GHz and Bi-RF, respectively. 
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In the right-hand side, the Bi state (i £ [1 , 20] ) is rep- 
resented by the product of the electron (mj = ±1/2) 
and nuclear spin (—9/2 < mj < 9/2) states with the 
total spin z-component of — 5 < M < 5 is represented 
on the basis of the electron (ms) and nuclear (mi) spin 
z-component eigenstates.- 2 ^. Note that there are two dif- 
ferent Bi eigenstates for one particular M, except for 
M = ±5. The mixing angle 4>i(Bo) depends on the Bi 
state and is a function of the parameters in the Hamil- 
tonian, Eq. (1), as explicitly described in Ref. [2(j. De- 
pending on the state of PRC (fi = ±1/2, denoted by 
fi =t, J.), each term in Eq. (5) gives contribution to the 
recombination rate in terms of spin parallel (i? p ) or anti- 
parallel (Rap) pair: 

Ri A (B ) = R p cos 2 &(B ) + R ap sin 2 fo(B ), (6) 

Ri,i(Bo) = i? p sin 2 4>i(B Q ) + i? ap cos 2 (7) 

Then, only the recombination associated with the pure 
states, i = 10 (M = —5) and i = 20 (M = +5), have sin- 
gle components that are strictly parallel or anti-parallel, 
while the other states have a mixture of the two com- 
ponents. This, in combination with Eq. (4), is the rea- 
son why the highest- and lowest-field lines, which involve 
the Bi state 1 10) or |20), are stronger than the other 
lines at X-band resonance and exclusively strongest at 
the radio frequency resonance. We used this model to 
perform the simulation of the X-band spectrum in Fig. 
1(b) and the 200 MHz spectrum in Fig. 3(a). The ratio 
Rp/Rap = 0.01 has led to good agreement with the exper- 
iments and is comparable to the recently reported value 
15 /is / 2 ms = 0.0075 for the phosphorus donor in Sii^S 
Note that, for the X-band spectrum we also take into 
account that an EPR transition line that is "forbidden" 
in the high- field limit overlaps with each EPR- "allowed" 
one except for the highest- and lowest-field lines that in- 
volve one pure stated 

Finally, we shall demonstrate the tunability of such 
"pure-state" transitions to the energy comparable to the 



superconducting qubits. As shown in Fig. 4 (a), an addi- 
tional 8.141 GHz microwave excitation in the same SDR 
method allows for successful excitation and detection of 
the EPR transition between Bi |1) and 1 20) levels at low 
field (Bo = 30 mT). Although it is preferred to achieve 
Bq < 10 mT for the coupling with the superconducting 
qubit, Fig. 4 shows clearly the flexibility to tune the 
energy difference between up and down states of the Bi 
electron spin. It is also possible to tune the supercon- 
ducting flux qubit to match the energy between \R) and 
\L) states with the Bi transition frequency separating |1) 
and 1 20) states. The coupling strength between the flux 
qubit and Bi is expected in the range of 1 — 100 kHz. 4 
Figure 4(b) shows a result of a similar experiment but 
with two excitation frequencies generated by two coils 
perpendicular to each other. The second coil was used to 
irradiate 100-MHz excitation frequency. The Bi-RF and 
Bi-8 GHz resonance lines are observed together. Hence, 
this experimental setup allows us to perform SDR-MR 
with two arbitrary excitation frequencies. 

In summary, we have obtained the electron paramag- 
netic resonance spectra of a small number (5 x 10 11 ) of Bi 
donors in Si at low magnetic field (6 — 110 mT). The de- 
tection was based on the measurement of the sample pho- 
toconductivity which changed significantly at the time of 
resonance due to specific spin-dependent-recombination 
phenomena. The spin-dependent-recombination process 
takes place via coupling of the electron spins between Bi 
donors and nearby paramagnetic recombination centers. 
The relative intensity of each resonance line has been 
described well by a spin-dependent-recombination model 
based on the mixing of Bi donor electron and nuclear 
spins. 
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Scientific Research and Project for Developing Innova- 
tion Systems by MEXT, FIRST, and JST-EPSRC/SIC 
(EP/H025952/1). 
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